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Fifteen months of cumulative CoGeNT data are examined for indications of an annual modula- 
tion, a predicted signature of Weakly Interacting Massive Particle (WIMP) interactions. Presently 
available data support the presence of a modulated component of unknown origin, with param- 
eters prima facie compatible with a galactic halo composed of light-mass WIMPs. Unoptimized 
estimators yield a statistical significance for a modulation of ~2.8cr, limited by the short exposure. 
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CoGeNT employs P-type Point Contact (PPC) germa- 
nium detectors [IrlSi to explore light-mass WIMP dark 
matter models. By virtue of their low electronic noise, 
PPCs are ideal for searches in the mass range m x < 10 
GeV/c 2 . Their ability to reject background events taking 
place on detector surfaces is described in jlj. Prompted 
by the interruption in data-taking imposed by a recent 
fire in the access shaft to the Soudan Underground Lab- 
oratory (SUL), we have examined existing CoGeNT data 
for the presence of an annual modulation in the WIMP 
interaction rate The characteristics of this detector, 
data analysis and background interpretations are treated 
in a recent Letter [lj, to which the reader is referred. 

Underground installation of this PPC at SUL took 
place on August 21, 2009. Following an upgrade to the 
data acquisition to allow discrimination of surface events, 
data-taking started on December 4th, a date close to the 
minimum in rate expected from the annual modulation 
effect Data-taking was interrupted exclusively over 
the periods of February 9-15, March 15-20, and Octo- 
ber 5-7 of 2010, for inspection of a higher-energy region 
of the spectrum and two planned general power outages 
at SUL. The dataset presented here ends on March 6 of 
2011, spanning 458 days, of which 442 were live. On this 
date the detector was stopped for another planned out- 
age. Outage-related computer problems delayed normal 
operation until two days before the fire (March 17, 2011). 

Surface background- and microphonic-rejection cuts on 
these data are as in [lj], and constant in time. An alter- 
native analysis [f| results into very similar conclusions, 
to be treated elsewhere Q. In the following discussion, 
the astrophysical halo parameters in 0, H| are used. 



Fig. 1 (top) displays spectral peaks appearing at the 
K-shell binding energy of daughters of cosmogenically- 
activated radioisotopes decaying via electron capture 
(EC) [lj]. The bottom panel in the same figure zooms 
in to the lower energy region of the spectrum, down 
to the ^0.5 keV ee (keV electron equivalent or ioniza- 
tion energy) threshold. A fit to the evolution of the 
K-shell peaks returns excellent agreement with the ex- 
pected half-lives and allows a prediction for the initial 
abundance of these radioisotopes with individual uncer- 
tainty of O(10)%. The ratio between L-shell and K-shell 
EC is well-defined for these isotopes, both theoretically 
and experimentally [t| . It can be used to generate a pre- 
diction, devoid of any free parameters, for the intensity 
of the L-shell peaks in this lower-energy spectrum (Gaus- 
sians in Fig. 1, bottom). After subtraction of this pre- 
dicted L-shell EC contribution and of a constant spectral 
component, and following the correction for the combined 
trigger and software-cut efficiency [l| , an exponential-like 
irreducible background of events taking place in the bulk 
of the crystal is observed (Fig. 1, inset). These events 
are individually distinct from any identifiable source of 
noise, display the characteristics (rise and decay time) 
of radiation-induced pulses, and have a uniform diur- 
nal distribution. Predicted signals from example com- 
binations of light- WIMP masses and spin- independent 
WIMP-nuclei couplings are also shown as a reference. 

Black dots in the inset of Fig. 1 represent the irre- 
ducible spectrum obtained by stripping of the L-shell 
predictions and the flat background level in the region 
~2-4.5 keVee- Unfilled circles are obtained when a free 
overall normalization factor multiplies the envelope of 
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FIG. 1: Top: Uncorrected (i.e., prior to threshold efficiency 
correction) spectrum displaying all expected K-shell EC cos- 
mogenic peak positions. The dotted histogram shows the 
spectrum before rejection of surface background events. Bot- 
tom: Uncorrected low-energy spectrum following removal of 
surface events. Dotted Gaussian peaks show the predicted 
L-shell EC contribution, devoid of any free parameters (see 
text). A dashed line traces their envelope. A second dashed 
line indicates the combined threshold efficiency (trigger + 
software cuts) [lj, an arrow pointing from it to the right scale. 
Inset: Spectra corrected by this efficiency and stripped of L- 
shell contribution and flat background component. Examples 
of light WIMP signals are overlapped on it (see text). 

the individual L-shell predictions in a background model 
containing this envelope, an exponential and a constant 
background. The resulting best-fit indicates a L-shell 
contribution just 10% short of the nominal prediction, 
well within its uncertainty. Fig. 2 shows the region of in- 
terest (ROI) obtained when these irreducible spectra are 
fitted by a sample model containing signals from WIMPs 
of mass m x and spin-independent coupling a si, and a 
free exponential background. As in [![, this ROI is de- 
fined by the upper and lower 90% C.L. intervals for the 
best-fit u si, whenever the lower interval is incompatible 
with a null value. This ROI is meant to direct the eye 
to the region of parameter space where the hypothesis of 
a WIMP signal dominating the irreducible background 
events fares best, but it does not include astrophysical or 
other uncertainties listed next. Reasonable uncertainties 
in the germanium quenching factor employed (Fig. 4 in 
0, can shift this ROI by - ±1 GeV/c 2 . The present 
uncertainty in the fiducial bulk volume of this detector 
is o(io)% Q. Departures from the assumption of a con- 
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FIG. 2: ROI extracted from the irreducible spectra in Fig. 1 
(inset) under consideration of a light-WIMP hypothesis. A 
small dotted line bisects it, approximately separating the do- 
mains favored by the black dot (left) or unfilled circle (right) 
spectra in Fig. 1. ROI definition and uncertainties able to 
shift it are described in the text. The DAM A/LIBRA ROI 
includes present uncertainties in its position [llj, with the 
exception of ion channeling [T3 ]. conservatively assumed to 
be absent. Solid and dotted lines are CDMS limits from [llj 
and 0], respectively. A dashed line corresponds to recent 
XENON100 claims [a|. U ncertainties in these constraints and 
those by XENON10 Tpl] are examined in [13, [3. 

stant background in the model above can also displace 
this region. A modest contamination of the spectrum by 
surface events next to threshold 0, @ would shift this 
ROI to slightly higher values of m x and lower asi- The 
additional exposure collected since |l[ results in a much 
reduced CoGeNT ROI, one in the immediate vicinity of 
the parameter space compatible with the annual modu- 
lation effect observed by DAMA/LIBRA QHG3. This 
region of a si, m x space is populated by the predictions 
of several particle phenomenologies. The reader is di- 
rected to references in jl[ and recent literature for ex- 
amples. The same region has received recent attention 
within the context of dark matter annihilation signatures 
at the center of our galaxy, and anomalies in accelerator 
experiments Fig. 2 also displays limits from other 

searches, a subject treated again below. 

A search for a WIMP-induced annual modulation in 
dark matter detector data requires an exceptional low- 
energy stability in the device. Fig. 3 shows that these 
conditions are present for CoGeNT. The top panel dis- 
plays daily averages in the detector electronic noise. Ex- 
cessive excursions in this parameter would affect the sta- 
bility of the detector threshold. These are not observed. 
Precautions are taken to ensure that this noise is as sta- 
ble as possible: for instance, by automatically refilling 
the detector liquid nitrogen Dewar every 48h, the crystal 
temperature and its associated leakage current are held 
as constant as possible. The second panel shows the sta- 
bility of the trigger threshold, derived from the difference 
between the daily average baseline DC level in the trig- 
gering channel and a constant (digitally fixed) discrimi- 
nator level. The small excursions observed correspond to 
a temperature drift in the digitizers (NI 5102) and shap- 
ing amplifier (Ortec 672) of ~ 1°C. These small instabil- 
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ities do not result in any minor smearing of the energy 
resolution, given that the amplitude of each event is refer- 
enced to its individual pre-trigger DC level. The effect of 
this small baseline drift should instead be envisioned as a 
maximum shift of the threshold efficiency curve in Fig. 1 
by about ±10 eV. The third panel shows the calculation 
of by how much such a shift can affect the counting rate 
in the region 0.5-0.9 keV ee . This calculation includes the 
exponential spectral shape observed there. The correc- 
tion is referenced to the date of the threshold efficiency 
calibration employed (small arrow in Fig. 3) and found to 
be negligible at less than 0.1%. This correction would be 
larger for events below 0.5 keV ee , not considered here, 
and even smaller for count rates in broader energy re- 
gions starting at 0.5 keV ee . The fourth panel indicates 
the magnitude of the correction required to account for 
the exponential decay of L-shell EC radioisotopes, prior 
to an annual modulation analysis. This correction affects 
the 0.4-1.6 keV ee region (Fig. 1), where a light- WIMP can 
express a modulated signature. The uncertainties in this 
correction, indicated in Fig. 3 in parentheses, are mod- 
est even at the present exposure. A direct comparison of 
these predictions with the low-energy spectrum, as done 
in Fig. 1, demonstrates that they are robust. 

Fig. 4 shows the temporal rate behavior in several spec- 
tral regions following the correction for L-shell EC activ- 
ity, when applicable. Even with the present limited ex- 
posure, a noticeable annual modulation can be observed 
in the energy region encompassing the exponential rise 
in irreducible background in Fig. 1. The statistical sig- 
nificance for this modulation is presently maximal for an 
energy bin ranging from threshold (0.5 keV ee ) to an up- 
per bracket anywhere in the interval ~2. 0-3.0 keV ee (the 
amplitude of a WIMP-induced modulation is expected 
to be maximal towards its spectral endpoint [HI, l2fjj|). 
No indication of a modulation is observed above this en- 
ergy, nor for events rejected as surface backgrounds in the 
same spectral region where it is the largest for bulk events 
(bottom panel in Fig. 4). The "internal clock" provided 
by the most intense K-shell peak at 10.37 keV ee exhibits 
no modulated deviation from an exponential decay, nor 
any significant time-dependent changes in mean energy. 

Ad hoc methods able to extricate a modulated com- 
ponent with maximal sensitivity [20| have not been at- 
tempted. A few considerably less sensitive indicators of 
modulation significance are presently offered. For in- 
stance, for the region 0.5-3.0 keV ee and (unoptimized) 
binning depicted in Fig. 4, a straightforward analysis re- 
veals a reduced chi-square x 2 /d.o.f.=7.8/12 (80% C.L. 
for acceptance) for the best-fit modulation. The null 
hypothesis (absence of a modulation) fares considerably 
worse at x 2 /d-o.f.=20.3/15 (84% C.L. for rejection). The 
likelihood ratio test indicates that the modulation hy- 
pothesis is preferred over the null hypothesis at 99.4% 
C.L. (2.8er). Intriguingly, the best-fit values for the three 
modulation parameters (16.6±3.8% modulation ampli- 
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FIG. 3: Assessment of the stability of the CoGeNT PPC at 
Soudan (see text). First panel: daily average in detector elec- 
tronic noise (shaping time 10 [is). Second panel: stability of 
the trigger threshold level. Third panel: negligible correction 
to the counting rate in the region 0.5-0.9 keV ee induced by 
it. Fourth panel: expected counting rate in this same region 
originating in L-shell EC. The observed stability augurs well 
for WIMP modulation searches using next-generation PPCs 
like those planned for the upcoming expansion of CoGeNT 
(C-4), MAJORANA, GERDA, and CDEX. 



tude, period 347±29 days, minimum in Oct. 16±12d) 
do not fall far from the predictions provided by the 
WIMP hypothesis (a calculable amplitude, a yearly pe- 
riod and minimum amplitude in late Nov. to early Dec. 
21 1 ) . The most uncertain of the se p redictions [2lT |22| , 
the amplitude, is derived as in 119], including the de- 
pendence on astrophysical parameters, target, thresh- 
old, WIMP mass, etc. For m x =7 GcV/c 2 an expected 
value of 12.8% is obtained (Fig. 4). While the apparent 
DAMA/LIBRA modulation is weaker (—2% of the low- 
energy rate), its enhancement for CoGeNT is expected in 
most light- WIMP scenarios 0, [23[ . If these predictions 
are accepted, the Monte Carlo probability of better si- 
multaneous agreement with them than that provided by 
the best-fit values above, is small at 0.7%, for simulated 
random fluctuations around the 0.5-3.0 keV ee mean rate. 

In addition to the basic estimators presented here, 
time-stamped CoGeNT data are available by request. 
These can be used not only for alternative analyses of 
significance, etc., but also to investigate non-cosmological 
effects that might generate this modulation and by ex- 
tension that observed by DAMA/LIBRA. In this respect, 
the muon flux at SUL varies seasonally by ±2%, and 
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FIG. 4: Rate vs. time in several energy regions (the last bin 
spans 8 days). A dotted line denotes the best-fit modulation. 
A solid line indicates a prediction for a 7 GeV/c 2 WIMP in 
a galactic halo with Maxwellian velocity distribution. Back- 
ground contamination and/or a non-Maxwellian halo can shift 
the amplitude of this nominal modulation (see text). Dotted 
and solid lines overlap for the bottom panels. 



radon levels by a factor ~4 [2J] . Muon-coincident events 
constitute a few percent of the low-energy spectrum [l| , 
limiting a muon- induced modulated amplitude to <<1% 
Q . Rejection of veto-coincident events does not alter the 
observed modulation. Radon displacement via pressur- 
ized LN boil-off gas is continuously maintained at 2 1/min 
within an aluminum shell encasing the lead shielding [25| . 
A radon-induced modulation would be expected to affect 
a much broader spectral region than observed [2fj| . 

The CDMS collaboration has recently claimed [3] to 
exclude a light- WIMP interpretation of CoGeNT and 
DAMA/LIBRA observations. Uncertainties affecting 
this claim are discussed in [13, H3- Observations from 
XENON10 and XENONIOOQ have been used to 
claim a similar rejection of light- WIMP scenarios. Un- 
certainties affecting these searches are examined in [l8| • 

In conclusion, presently available CoGeNT data favor 
the presence of an annual modulation in the low-energy 
spectral rate, for events taking place in the bulk of the 
detector only. While its origin is presently unknown, the 
spectral and temporal information are prima facie con- 



gruent when the WIMP hypothesis is examined: in par- 
ticular, the WIMP mass region most favored by a spectral 
analysis (Fig. 2) generates predictions for the modulated 
amplitude in agreement with observations, modulo the 
dependence of this assertion on the choice of astrophysi- 
cal parameters and halo velocity distribution j2lM23l l28| . 
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